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We study the ‘flux noise’ spectrum of random-bond quantum Heisenberg spin systems using a 
real-space renormalization group (RSRG) procedure that accounts for both the renormalization of 
the system Hamiltonian and of a generic probe that measures the noise. For spin chains, we And that 
the dynamical structure factor Sq{f), at finite wave-vector q, exhibits a power-law behavior both at 
high and low frequencies /, with exponents that are connected to one another and to an anomalous 
dynamical exponent through relations that differ at T = 0 and T = oo. The low-frequency power- 
law behavior of the structure factor is inherited by any generic probe with a finite band-width and 
is of the form l//“ with 0.5 < a < 1. An analytical calculation of the structure factor, assuming a 
limiting distribution of the RG flow parameters (spin size, length, bond strength) confirms numerical 
findings. More generally, we demonstrate that this form of the structure factor, at high temperatures, 
is a manifestation of anomalous diffusion which directly follows from a generalized spin-diffusion 
propagator. We also argue that l//-noise is intimately connected to many-body-localization at 
finite temperatures. In two dimensions, the RG procedure is less reliable; however, it becomes 
convergent for quasi-one-dimensional geometries where we And that one-dimensional l//“ behavior 
is recovered at low frequencies; the latter configurations are likely representative of paramagnetic 
spin networks that produce l//“ noise in SQUIDs. 


I. INTRODUCTION 

Disordered interacting spin systems display a variety 
of physical phenomena: the spin glass transition [1-5] 
and slow relaxation [6-9], many-body localization tran¬ 
sition [10-14] and a breakdown of ergodicity, Griffiths 
(rare-region) effects [15-18], strong-randomness fixed 
points [19-21], spin-liquid states [22] and even excited 
topological states [23]. A great deal of our understand¬ 
ing of these phenomena has come from the develop¬ 
ment of the real-space renormalization group (RSRG) 
method [24, 25]. Previous applications of the method fo¬ 
cussed on studying the low-temperature thermodynamic 
properties of disordered spin chains. More recently, the 
approach has been applied to study both high temper¬ 
ature and dynamical properties—this is based on the 
identification that the protocol, more generally, involves 
the creation of local integrals of motion, i.e., the high- 
energy modes eliminated during the RG process are 
pieces of approximate many-body eigenstates of the sys¬ 
tem [19, 26, 27]. Using such ideas, theoretical analy¬ 
sis of the low-temperature optical conductivity of Anti- 
Ferromagnetic spin chains of various kinds was carried 
out in Ref. [19] while that for Ising spin chains was per¬ 
formed computationally to identify a transition between 
various infinite-temperature many-body localized phases 
in Ref. [26]. 

The work presented here builds on and extends the 
RSRG program to address a crucial question in the dy¬ 
namics of disordered systems—how do they generate 
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scale-invariant l/a;“ (or l//“) noise? Numerous ex¬ 
periments [28-34] on Superconducting Quantum Inter¬ 
ference devices (SQUIDs) observe a flux noise with a 
spectrum N(u}) ^ l/a;“ whose magnitude is nearly tem¬ 
perature independent, but exponent a changes smoothly 
with temperature [29, 30]. The noise likely originates 
from fluctuating electronic spins localized on the metal- 
insulator (the conducting strip and the substrate) inter¬ 
face of the SQUIDs [28, 35, 36]. Such spins can inter¬ 
act via oscillatory RKKY (Heisenberg) interactions [36] 
whose sign flips on the order of the Fermi-wavelength 
and is effectively random at the scale of separation of 
the spins, i.e., interactions can be both Ferromagnetic 
(F) or Anti-Ferromagntic (AF). The dynamics of such 
spins, and how they generate the observed noise spec¬ 
trum is less understood. In particular, it was posited that 
the two-dimensional surface spins exhibit regular diffu¬ 
sion [36, 37] and that the l/w noise spectrum arises in a 
limited frequency range owing to specifics of the geom¬ 
etry of the probe coupling. Such explanations predict a 
large frequency lower bound to the 1 /lo form of the noise 
spectrum which is not observed in experiments. More 
importantly, it was assumed that the disorder is self¬ 
averaging (which leads to diffusion). A central finding 
of our work is that, for sufficiently strong initial disor¬ 
der, this is not the correct conclusion for one-dimensional 
and two-dimensional ‘strips’ of Heisenberg spins—we find 
that such interacting spin networks generically flow to¬ 
wards strong-randomness fixed points where they exhibit 
anomalous diffusion and, moreover, this entails a flux 
noise that intrinsically exhibits a 1 /cli“ (a < 1) noise 
spectrum. 

The RSRG protocol we develop allows us to directly 
compute the noise spectrum N{ui) measured by a probe 
that couples (through arbitrary geometrical factors) to 
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FIG. 1. (a) Phase diagram of the strongly disordered Heisen¬ 
berg spin chain (with 1/| J| distribution of initial couplings in 
a range |J| G ) extrapolated to finite temperatures us¬ 

ing RSRG results at zero and infinite temperatures; the color 
scheme encodes the variation of the dynamical exponent 2 as 
a function of the bias rji [initial proportion of F/AF (1/-1) 
bonds in chain], (b) The low-frequency noise power law a 
and (d) the dynamical critical exponent 2 are plotted against 
T^i at T = 0 (data-points and error bars) and T = 00 (flat 
line with dashed lines indicating error), for disorder strength 
D — 3 (see main text). At T = 0, the purely F system har¬ 
bors heavily damped spin wave excitations whose localization 
length diverges as 1/tu in the zero-frequency limit. The purely 
AF system, which is known to flow to a infinite-randomness 
fixed point, exhibits 1/u) noise. Intermediate biases (non- 
hatched region) flow to strong-randomness fixed points ac¬ 
companied by l/a;“ noise spectra. Hatched region is inacces¬ 
sible to RSRG. At high temperatures (T 3> a;), the l/a;“ noise 
spectra is generic to strongly disordered Heisenberg chains 
and can be interpreted in the context of generalized diffusion. 
The crossover from zero-temperature to infinite-temperature 
behavior occurs at T/oj ~ 2^°, where zq is the dynamical 
exponent found for the zero temperature system. 

the flux generated by spins interacting via disordered ex¬ 
change couplings, at both high and low temperatures. 
Our approach is elaborated upon in Sec. II. A numeri¬ 
cal application of the protocol suggests that a (strongly) 
disordered Heisenberg spin chain flows towards a strong¬ 
randomness fixed point characterized by an anomalous 
dynamical exponent 2 = 1//3 1,2. Moreover, this dy¬ 

namical exponent determines the form of the noise; a har¬ 
monic probe with wave-vector q measures the noise Sq{u}) 
(the dynamical structure factor) which shows a piece- 
wise power-law frequency dependence: Sq(uj) ~ 
for w < and Sq{uj) ~ l/a;“' for w > The 

precise values of these power laws are independent of 
the initial composition (proportion of F/AF bonds) of 


the spin chain if the distribution of couplings is suffi¬ 
ciently non-singular, but varies for more singular initial 
distributions such as those relevant for experiments on 
SQUIDs (and that we primarily consider)—for such dis¬ 
tributions, the low-frequency noise exponent varies in the 
range 0.5 < a < 1 (Fig. 1). [Note that any probe of finite 
spatial bandwidth will inherit the low frequency exponent 
a in the structure factor 5 ^( 0 ;).] Our numerical results 
for spin chains are discussed in Sec. III. 

While exponents 2 , a, and a' can take a range of 
values, they are tied by the non-trivial relations a' = 
1 + 2(1 — a), a = 1 — Ijz at high temperatures, and 
a' = 1 + 3(1 — a), a = 1 — 1/22 at low temperatures. 
In Sec. IV, we show how these relations can be obtained 
from a scaling form of the probability distribution func¬ 
tion Pf(af) that governs the distribution of F(AF) bonds 
with a given coupling at energy scales where the RG pro¬ 
cedure has converged. These relations are then verified 
numerically by performing a scaling collapse of Sq{u}) for 
a wide range of wave-vectors and frequencies (Fig. 4) 
at both T = 0, and T = 00 ; the success of the scaling 
collapse gives further credence to the validity of the RG 
procedure. 

While the above discussion pertains specifically to the 
Heisenberg spin chain, we show, in Sec. V, that the form 
of the dynamical structure factor carries over generally, 
at high temperatures, independently of dimensionality, 
to systems exhibiting anomalous diffusion. We propose 
a generalized-diffusion ansatz for the spin propagator, 
Gq{u}) = l/[—iw -\- g^/^/(a;/g^/^)], that can account for 
anomalous diffusion and show that it directly reproduces 
the limiting g— and w—dependencies of the structure 
factor. Thus, we conclude that at high temperatures 
(T 3> w), l/w“ noise is generic to a system exhibiting 
anomalous diffusion (accompanied by an anomalous dy¬ 
namical exponent). The generalized-diffusion approach, 
however, fails to explain the structure factor at T = 0. 
We discuss how this is a consequence of the failure of 
linear response at T = 0—we show that the system ex¬ 
hibits a divergent static susceptibility for even finite-g 
perturbations. 

The l/a;“ low-frequency behavior of the structure fac¬ 
tor also has the consequence that all (harmonic) spin cor¬ 
relations decay to zero in the long time limit as 
Conversely, when a = 1 precisely, even finite-g modes, 
which are not conserved are unable to equilibrate, sig¬ 
naling the onset of the many-body localization (MBL) 
transition. Thus, l/w noise is connected to the MBL 
transition. A hint of this is observed in our simulations at 
T = 00 , where, upon increasing the disorder strength, the 
low-frequency noise exponent approaches 1, but does not 
reach it. At T = 0, the noise exponent 1 is obtained in 
the case of a purely AF chain, where the RG is expected 
to flow towards m/imte-randomness fixed point [24, 38]. 
The status of ergodicity in our model is less clear. The 
RG constructs a macroscopic number of integrals of mo¬ 
tion which suggests a lack of ergodicity. However, these 
integrals of motion are not exact since the RG flows to- 
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wards a strong-randomness fixed point, as opposed to 
infinite-randomness where it does become probabilisti¬ 
cally exact. We anticipate future studies should elucidate 
this question. (While MBL implies a lack of ergodicity, 
the converse is not true; see for instance Ref. [39].) 

In Sec. VI we discuss the application of the RSRG 
procedure to two-dimensional spin networks. For square 
geometries, at high temperatures, the structure factor 
shows scaling behavior (although in a limited range of 
frequencies) according to the generalized-diffusion ansatz 
for a dynamical exponent z ^ 2. Note that regular dif¬ 
fusion is also associated with z = 2; however, in the 
generalized-diffusion framework, this implies that the 
finite-g structure factor has a 1 /^/ui low-frequency tail 
instead of the fiat spectrum expected in the case of reg¬ 
ular diffusion. We note that for square geometries, the 
accuracy of the decimation procedure does not improve 
significantly over the course of the RG. Therefore, we 
cannot rule out (or confirm) regular diffusion. In con¬ 
trast, for an elongated two-dimensional ‘strip’, we find 
that at the lowest frequencies, this l/Vw tail transitions 
into an anomalous l/a;“ tail (0.5 < a < 1) as in the 
one-dimensional case; thus, as the clusters become wider 
than the width of the two-dimensional strip, their collec¬ 
tive dynamics mimics that of a one-dimensional array of 
disordered spins. We end by summarizing our results and 
their connection to experiments in Sec. VII and propos¬ 
ing further possible extensions of our RSRG protocol to 
systems with anisotropic couplings. 


II. RSRG APPROACH FOR COMPUTING 
NOISE 

The model we study is the disordered Heisenberg spin- 
1/2 system with Hamiltonian H = Jij’^ 

are picked independently from a distribution PoiJ) with 
support over both negative and positive values of J. The 
probe measures a magnetic flux M = '^^giSi, and the 
noise measured by such a probe has a spectral form 
N{ijj) = J^{^[M(t)., M(0)] + ^}, that is, the Fourier trans¬ 
form of the autocorrelation of M{t). (Such an auto¬ 
correlation directly yields the x—, y—, z—component-wise 
noise auto-correlations due to spin-rotation symmetry of 
the model. Also note that the methodology presented 
here can be straightforwardly extended to the case where 
giS are vectors measuring a certain the projection of the 
local spin, which can be non-uniform. We do not treat 
this general case here for the sake of simplicity.) The 
RSRG approach and the evaluation of the noise spec¬ 
trum can be summarized succinctly as follows. We first 
find the 2 spins Sa and Sb that are most strongly cou¬ 
pled to one another; these spins tend to process rapidly 
around their total moment = Sa + Sb, which on the 
other hand, moves slowly if the constituent spins are rel¬ 


atively weakly coupled to the external neighbors. The 
total spin S^ is then kept as an effective spin with an ef¬ 
fective interaction with the probe and its neighbors, while 
the rapid precession of the constituent spins is counted 
towards the ‘noise’ at the frequency at which they pro¬ 
cess about 5”+. This procedure is repeated ad-infinitum. 
The approach is illustrated in Figs. 2 (a) and (b). 

It must be noted that the quantum-mechnical deriva¬ 
tion of the ground state RSRG rules for bond decima¬ 
tion in mixed AF/FM case for was first carried out in 
Ref. [40]. However, as mentioned above, to calculate the 
noise spectrum at arbitrary temperature, we need to ex¬ 
tend these rules to arbitrary composite states and pro¬ 
vide renomalization rules for the probe function gi. In 
order to keep the discussion self-contained and shed light 
on the physics of the problem, we discuss the derivation 
of these rules using a straightforward semi-classical ap¬ 
proach, that yields the quantum-mechnical results upon 
spin re-quantization. 


A. Derivation of RG rules 


To obtain the general bond decimation rules, let us 
first consider a set of two strongly interacting spins Sa 
and Sb with mutual coupling J and couplings Ja and 
Jb to external spins, such that J ^ Ja,Jb- Due to 
their large mutual coupling J, spins A and B precess 
about their combined moment S+ = Sa + Sb with a 


large frequency Aq = J 


The neighboring spins 


couple only to the slow motion of spins Sa and Sb 
which is the projection of these spins on to the slow com¬ 
bined moment S'+. Thus, the effective couplings Ja{b) 


54 


The 


are modified to Ja(^b) = Ja{b){Sa{b) ■ S+)/ 
quantum-mechanical result is obtained by re-quantizing 
the spins, i.e., |5 aP = sa(s41 + 1), |5bP = sb{sb + 1), 
and |5+P = s(s -I- 1). 

The quantum-mechanical interpretation of the above 
procedure is that, at every step of the RG, we solve the 
approximate Hamiltonian Hab = JSa-Sb + ... to first 
order in perturbation theory by projecting the system 


to a particular eigenstate of S+ with quantum number 
s. This effective spin has couplings and Jg deter¬ 
mined above with its neighbors. The frequency Aq with 
which the spins precess about the total moment S^- is 
split into 2 in the quantum-mechanical situation, and 
is given by the energy difference between total angu¬ 
lar momentum states s and s ± 1 of the coupled spins. 
The perturbation theory is controlled by the parameter 
Aa(_b)/Ao where Aa,^b correspond to a similar gaps 
generated by the combination of spins A and B with their 
neighbors. Note that in the case of singlet-formation 
(s = 0), second-order perturbation theory must be per¬ 
formed, eliminating spins A and B altogether and yield¬ 
ing an effective coupling between their neighboring spins, 
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FIG. 2. (a) Two-spin dynamics. Two spins A and B strongly coupled to one another revolve around their net moment 

5+ = Sa + Sb with a large frequency J|5'+|, while S+ moves slowly owing to weak interactions J'a and J^. S_ = Sa — Sb 
is effectively decoupled from the slow dynamics of 5+ and results in the ‘noise’ evaluated at this RG step, (b) Demonstration 
of RG rules. The strongly coupled spin-pair are decimated to form an effective spin (with quantized net moment s+), and 
Heisenberg couplings are renormalized from Ja(b) to Ja'(b')- The effective spin couples to the probe with amplitude <?+. If a 
singlet is formed, both spins are integrated out, providing an effective interaction J' between the neighboring spins, (c) Scaling 
behavior at T = oo (at T = 0, the bonds, composed of long-range singlets scale as the cluster size). Size of clusters n, their 
spin size s and the time-scale of their dynamics 1/uj are connected through the anomalous dynamical exponent 1//3. Harmonic 
probes (illustrated in blue) measure a low-(high-)frequency noise-exponent a or a' in the regimes where the probe length 1/q 
is smaller (greater) than the clusters whose dynamics they probe. 


J' = 2sa{sa + ^)JaJb/'^J {sa = Sb) [40]. Note that in 
higher dimensions, when external spin A (and/or B) may 
couple to both spins 1 and 2 undergoing a singlet forma¬ 
tion with couplings Ja,i and Ja .,2 respectively, the cou¬ 
pling Ja in this result should be replaced by Ja,i — Ja, 2 - 
The renormalization of the probe coupling is also im¬ 
mediate: if the probe couples with strength qa and gs to 
spins A and B, it couples to the effective spin with a 


strength g+ = gA^A ■ S+)/ 

{gA + gB)l‘2- + [gA — gB)/‘2‘ 


Sa 


Sa 


+ gB{SB-S+)/ 


Sa 


Sb 


/ 


^4 


Note 


that in the case of singlet formation, the spins are elim¬ 
inated altogether and the probe does not couple to the 
singlet. This will result in important differences in the 
behavior of the system at T = 0 (where singlets are pre¬ 
ferred) and T = oo (where singlets are unfavorable for 
entropic reasons). 


B. Evaluation of Noise at every RG step 


We want to evaluate the dynamics of the object 
M{t) = J2igiSi- As mentioned above, we do this in 
a step-by-step RG fashion. In particular, the flux can 
be partitioned into a slow part Ms{t) = b gi^i + 


g+S^ whose dynamics is determined in subsequent RG 
steps, and the remaining fast part Mp(t) = {gA — 


gB)!^ 


s. - Spis. • Sp)! 


, where S- = Sa — Sb 


and the combination in square brackets is the compo¬ 


nent of S- orthogonal to Sp. To zeroth order in pertur¬ 
bation theory (in Aa(b )/the noise spectrum N{uj) 
receives a contribution B{{[MF{t), Mf{0)]+)} from this 
step of the RG, where the brackets () correspond to a 
quantum mechanical expectation in an eigenstate |s,m) 
with |5'+P = s(s -|- 1) and an arbitrary projection quan¬ 
tum number m (there is no preferred axis for the total 
moment 5'+). 


Note that, S- produces transitions from states 
|s, m) to states |s±1,to±1). Due to the full ro¬ 
tational symmetry of the problem, the transition 
rates are independent of m, and the transition 
frequencies only depend on s. Hence, the noise 
associated with these contributions is N{uj) = 
{gA - gB)'^/^[M{s,t)S{uj-uj^)+M{s,i)S{u} + u}i)], 

where is the frequency of the transition from the 

angular momentum state s to s -I- 1 (s to s — 1) and 
M(s,t (/)) is the accompanying matrix element of this 
transition. In what follows, we will refer to the factor 
{gA — gB)"^ 1^ as the ‘probe form factor’ since this factor 
depends on the precise details of the probe, and the 
matrix elements M(s,t (i)) corresponding to the spin 
operators will be referred to as the ‘noise amplitude’. 
The precise evaluation of the noise amplitude is discussed 
in Appendix A. 







































5 


III. NUMERICAL SIMULATIONS IN ONE 
DIMENSION 

A. Implementation 

As a first check of the efficacy of the RG rules, we per¬ 
form a direct comparison of the complete energy spec¬ 
trum determined by exact-diagonalization and the RG 
procedure for small random spin chains. A typical run 
is shown is Fig. 3; the RG result appears to be in good 
agreement with the exact-diagonaliation result. (More 
quantitative checks of the RG procedure are discussed 
in Appendix E.) For the determination of the noise and 
the structure factor, we perform numerical simulations 
for spin chains of size L = 15000 at two fixed temper¬ 
atures T = 0 and T = oo. The majority of our sim¬ 
ulations assume distribution of magnitude of the initial 
couplings Po(|<.^|) is of the form Pod-.^!) ~ with val¬ 
ues in the range |J| G This choice is motivated 

by the observation that a system of spins scattered ran¬ 
domly in d dimensions interacting with a 1/r^ RKKY 
interaction, has a distribution Po(| J|) ~ 1/| the 

additional factor of d/S is not systematically considered 
since it is found to not affect the results qualitatively. (In 
Ref. [40], Westerberg et al. work with an array of initial 
distributions and conclude that distributions more sin¬ 
gular than Pod J|) ^ 1/1 J|°-^ flow to non-universal fixed 
points at T = 0.) Another simplification we make is to 
consider interactions only between nearest-neighbor spins 
even though spins interact with all other spins via the 
RKKY mechanism; this is well justified because l/r^ in¬ 
teractions are sufficiently short-ranged in one and quasi- 
one dimensional cases we primarily consider. We also 
perform simulations with a uniform distribution (range 
I J| G [0,1]). We choose a finite proportion of these cou¬ 
plings to be AF(F) and characterize this initial ‘bias’ by a 
variable G (—1,1): r]i = -|-1(—1) corresponds to a sys¬ 
tem with purely AF(F) couplings. To perform the zero 
temperature calculations, at every RG step, we choose 
the effective spin to reside in a state with s = |s .4 — Sb| 
( s = sa + SB ) for spins A and B that are coupled by 
a AF(F) bond. At infinite temperature, the quantum 
number s is chosen probabilistically according to the de¬ 
generacy 2 s + 1 associated with the state. 


The choice of the probe is another free parameter to 
be considered in our problem. The most natural choice 
of probe is gi = cos{qi), that is, a harmonic probe that 
measures the dynamical structure factor Sq{uj) at a given 
wave-vector q. A generic probe simply measures a noise 
= J2q \9q\^Sq{u}). 



FIG. 3. RG-determined spectrum (red) compared with exact 
digaonlization determined spectrum (blue) for a particular 12 - 
site random spin chain with initial 1 /| J| distribution, disorder 
strength D = 3, and equal mix of F/AF bonds. 


B. Convergence 

It was pointed out in Ref. [40] that the uniform distri¬ 
bution flows to a universal fixed point (at zero tempera¬ 
ture) with a final bias rif « 0.26 irrespective of its initial 
bias rji^ while more singular distributions, such as the 
ones we primarily consider do not flow to that particular 
fixed point. While we recover this result for initially uni¬ 
form distributions in our simulations, we note that for the 
1/|J| distributions we consider, the system also flows to 
fixed points (with a stable final bias 77 / G (0,0.3), see Ap¬ 
pendix D) but one that varies slightly depending on the 
value of the initial bias rji, disorder strength D, and tem¬ 
perature. Regardless of the precise values of the initial 
bias, the form of the distribution, or the temperature, we 
find that the effective gaps (and the corresponding effec¬ 
tive couplings, as well) flow to a power law distribution 
Pa(A) ~ 1/A''', 7 < I. This validates the RG proce¬ 
dure because a singular power law distribution results in 
a typical value of the gap A that is significantly smaller 
(by a factor e^Ai- 7 )) than the maximum gap Aq and 
guarantees a separation in energy scales between nearby 
regions. We also note that while the RG protocol re¬ 
sults in larger and larger spin sizes over the course of the 
RG, these spins also couple more weakly to their neigh¬ 
bors (see Appendix H) and this ensures that the RG can 
converge to a fixed point with stable scaling properties. 

The final bias 77 / and 7 display the following trends: 
i) at T = 0, for a uniform distribution, we find in con¬ 
firmation with Ref. [40], that rjf and 7 are independent 
of the initial bias; ii) at T = 0 , for 1 /|J| initial distri¬ 
bution, these quantities depend on the initial bias. In 
particular, 7 goes from 1 to 0 as the initial bias is varied 
from a purely AF to purely F chain, signaling the break¬ 
down of the RG procedure for the purely F chain; and 
iii) at T = 00 , for both forms of distributions, 77 / « 0 
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FIG. 4. Scaled dynamical structure factor (a) Sg(oj) = at T = 0, (b) S^{ui) = x ~^^^at T = oo 

for q G 27r/[40(yellow), 80(cyan), 120(green), 400(blue), 800(red)], and x = (27r/40)/5; the results shown correspond to a system 
of an equal proportion of F/AF bonds, and disorder strength D — 3; (c) scaling relations at T = 0 and T = oo between the 
power laws a and a' determining the high- and low-frequency dependence of Sq{oj) as given in Eq. (1). a and a' were plotted 
for D = 3 and r/i = [1, 0.75,0.5,0.25, 0,—0.25,—0.5,—0.75] for T = 0, while, at T = oo, the results are independent of rji, and 
the data-points correspond to D = [1,3,9]. Decreasing D and decreasing rji correspond to smaller values of a. At T = oo, 
a' K, 1 -\- 2(1 — a), while at T = 0, a' ~ 1 -I- 3(1 — a); these relations are plotted as dashed lines. The optimal scaling collapse 
in (a) (at T = 0) and (b) (at T = oo) was found for P = 0.40 and /3 = 0.36 respectively. These values are in accordance with 
analytically determined relation to the low-frequency noise exponent a, i.e., p = 2(1 — a) at T = 0 and (1 = 1 — a at T = oo. 


always, and 7 is independent of the initial bias but de¬ 
pends weakly on disorder strength D for the 1/| J| initial 
distribution. Further details on these observations can 
be found in Appendix E. 

[Additional validity checks : i) simulations for classical 
spin chains using classical RSRG rules at T = 00 were 
found to compare well with the quantum T = 00 RSRG 
results, see Appendix L ; ii) the structure factor obtained 
from the classical RG results was found to be in good 
agreement with that obtained from direct integration of 
Bloch-dynamics of the classical spins, see Appendix Mj. 


C. Form of the structure factor 


A most interesting facet of the simulations is the 
emergence of a finite, anomalous dynamical exponent 
z = 1/P ^ 1,2. This dynamical exponent dictates that 
the size of clusters, n, and the time-scale of their dynam¬ 
ics, uj~^, scale as n ~ uj~^ [as illustrated in Fig. 2 (c)]. 
Furthermore, this scaling gives rise to a piece-wise power- 
law behavior of the frequency dependence of the structure 
factor Sq{uj) for w and uj ^ corresponding 

to whether the probe’s period of oscillations 27r/g is much 
greater than the cluster size or much smaller. Moreover, 
the exponent /3 is directly connected to the exponent a 
in the noise spectrum N{uj) ^ l/w“. 

The full form of the dynamical structure factor Sg{uj), 
defined as Sq{uj) at T = 0 and S^{uj) at T = 00 can be 
extracted using a scaling collapse [Fig. 4 (a), (b), see also 
Appendix C for verification of g—dependent scaling] and 
is found to be 




I l/a;i-/5/2 < qVP, 

|^^2/^l+3/3/2 

( 1 ) 


where go(x) and gooix) are arbitrary scaling functions 
whose limits a; <C 1 , a; ^ 1 have been verified numeri¬ 
cally. (The limiting forms of the structure factor are also 
derived using an analytical approach in Sec. IV.) The 
low-(high-)frequency power-law exponent in both cases is 
identified as a {a'). As we had mentioned before, Sq{ui) 
shows a power-law frequency tail of the form l/a;“. It 
is curious to note that, at high temperatures, the noise 
exponent deviates from I by the inverse of the dynamical 
exponent, i.e., a = 1 — Ijz while at low temperatures, 
the appropriate relation isa = l — l/2z. The frequency 
dependence of the noise at higher frequencies (character¬ 
ized by exponent a') is also related in a non-trivial way 
to the dynamical exponent that differs at low and high 
temperatures. 

Using the scaling form of the structure factor, one can 
easily find that a typical probe that measures the noise 
N{oj) = J2q \ 9 q\‘^Sq{u}) will not observe any lower bound 
in the l/w-like behavior of the noise. In particular, for 
frequencies w ^ qq^^, where qg characterizes the probe’s 
resolution so that gq{q < go) ~ const., the noise spec¬ 
trum will inherit the low-frequency behavior of this finite- 
g structure factor. 

As for the precise values of these exponents, we find 
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at T = 0, for the 1/| J| initial distribution, that the low- 
frequency noise-exponent a varies from 1 (/3 = 0) for the 
AF chain to approximately 0.5 (/3 = 1) for the nearly F 
chain. The AF result of a « 1 can be understood from 
the point of view that the purely AF Heisenberg spin 
chain flows to an infinite-randomness fixed point with a 
distribution of couplings T’(T) ~ 1/J. In the nearly AF 
case, these couplings also directly give the gap A = J(l-|- 
|si ~ 521 ) ~ >7 since the system forms only singlets and 
the spins don’t grow in magnitude over the course of the 
RG. Since the noise generation happens at the frequency 
A (in particular, at every RG step, the maximum gap 
Ao), it simply inherits the power law of the coupling 
distribution. (The connection between the noise power a 
and the power law of the gap distribution is detailed in 
Appendix F) 

As we introduce more ferromagnetic bonds into the 
spin chain, the low-frequency noise exponent gradually 
approaches 0.5, and the distribution of effective gaps be¬ 
comes less singular (accompanied by a reduction in the 
fidelity of the RG). For an initial bias r\i ~ —0.75, the dis¬ 
tribution of gaps is almost uniform indicating a failure of 
the RG. Note that, at T = 0, a —>■ 0.5 is accompanied 
by a dynamical exponent z = 1//3 —>■ 1 (recall that ex¬ 
ponents are related at T = 0 via a = 1 — ll2z), and one 
can show (see Appendix J) that in such a case the RSRG 
fails due to the proliferation of faraway resonances. 

To gain more insight into the behavior of nearly F 
chain, we perform a Holstein-Primakoff (see Appendix K) 
expansion on top of the fully polarized ground state of 
a purely F chain. Unlike Eq. (1), the structure for the 
purely F spin chain is peaked at a finite frequency oj 
for wave-vector q that scales as a; ~ as expected for 
magnons. However, these magnons are heavily damped; 
the width of the magnon peak also scales as q^. Due to 
the large damping, these spin waves are ‘localized’ in the 
sense that they exhibit a finite inverse participation ratio 
[given by nf / Ui)^ ; Ui is the spin wave density at 
site i of the chain]. If we interpret this finite inverse par¬ 
ticipation ratio as a localization length, we find that it 
diverges in the low-frequency limit as ^(w) ^ l/w; this is 
in agreement with the observation that spin clusters scale 
as n{uj) ^ 1/uj (z = 1) for the nearly F chain in the RG 
at T = 0. (In two dimensions, the Holstein-Primakoff 
analysis shows that the spin-wave peaks are sharper and 
delocalized; the results are in agreement with effective- 
medium approximation based approaches [41, 42].) 

At infinite temperature, a is independent of the the 
initial bias, but depends weakly on the disorder D, ap¬ 
proaching 1 as disorder is increased [a crossing beyond 1 
is not observed, see Fig. 4 (c) where a and a' are plot¬ 
ted]. The range of the low-frequency noise-exponent is, 
as in the T = 0 case, 0.5 < a < 1. In Sec. V we will 
argue that the infinite-temperature results we obtain are 
natural for any system exhibiting anomalous diffusion. 
Moreover, we will see that the special case of a = 1 (or 
l/w noise, precisely) is singular enough that it leads to 
the absence of relaxation of non-conserved finite-g modes 


signaling many-body localization. It seems reasonable, in 
the light of these arguments, that the exponent a should 
approach 1 as the disorder strength D is increased [as we 
find, see Fig. 4 (c)]. 

We now recapitulate the findings for the strongly dis¬ 
ordered Heisenberg spin chain with initial coupling dis¬ 
tribution Pod'll) ^/\J\ il'^l S [l,e^]. At T = 0, 
we find that the system has three distinct regimes ac¬ 
cording to the initial bias: a) the purely AF spin chain 
which flows towards infinite-randomness and has a diver¬ 
gent dynamical exponent and l/ui noise spectrum; b) the 
purely F spin chain which has peaked spectral functions 
associated with heavily damped magnons whose size n(u}) 
diverges as l/w; and c) a wide range of mixed AF/F 
spin chains which flow towards non-universal strong¬ 
randomness fixed points (thus, accessible by RSRG) and 
whose structure factor is of the form in Eq. (1). These 
chains exhibit 1/uj°‘ noise for 0.5 < a < 1, which arises 
due to spin clusters whose size n{uj) diverges as 
with z > 1 ; this divergence is slower than 1/oj found for 
the purely E chain. At T = oo, the physics of the spin 
chain is independent of the proportion of AE/E bonds. It 
exhibits a l/a;“ noise spectrum, again with 0.5 < a < I 
where a increases as D is increased. Note that at T = oo, 
the dynamical exponent is related to a via the relation 
a = 1 — z, implying that z > 2 at T = oo. 

We can qualitatively glean the form of the structure 
factor at finite temperatures by considering a straight¬ 
forward extension of the RSRG protocol: we follow the 
zero-temperature RSRG rules for eliminating modes with 
frequencies w > T (since we expect these modes to be 
populated primarily in their ground state configurations), 
and infinite-temperature rules for eliminating modes at 
frequencies uj (since we expect these modes to be suf¬ 
ficiently excited). As we discuss below, such a protocol 
leads to the conclusion that the dynamics of the system 
resemble the T = oo (T = 0) behavior for frequencies 
w ^ T/ 2 ^“ (w 3> T/2^° ), where zg is the dynamical 
exponent of system (with a given initial bias rji) found 
at T — 0. Over the course of eliminating high-frequency 
modes (w > T), the gap distribution develops a power- 
law form with an exponent that is expected for the zero- 
temperature system. The zero-temperature RG proto¬ 
col continues until Ag « T. Beyond this scale, the RG 
continues with infinite-temperature rules, which begin to 
have an effect on the gap distribution only when a signif¬ 
icant fraction, say 1/2 of the spins have been eliminated. 
Using the dynamical exponent zg, the maximum gap Ag 
at this length scale is given by Ag = T/2^°. Thus, we ex¬ 
pect that the infinite-temperature form of the structure 
factor appears only for frequencies below w < T/2^°. 

Our findings are summarized in the phase diagram of 
Fig. 1. 


IV. SCALING APPROACH TO STRUCTURE 
FACTOR 

A. Scaling distribution Pf(af) 

We now explain the results in Eq. (1) by extending the 
arguments of Westerberg et al. (Ref. [40]) to evaluate 
the dynamical structure factor. They find that the RG 
flow eventually converges to a fixed point where the F/AF 
bias stabilizes and the probability distribution for F (AF) 
bonds is given by the scaling form Pp(^^)(A, Sl, Sr) ^ 

, SrA^^'^), where Aq corre- 

spends to the maximum gap at any step of the RG, while 
Sl, Sr are the left and right spins across the bond of fre¬ 
quency A. Note that, at T = 0, A is defined as the energy 
difference between first excited state from the ground 
state of the F/AF spin pair, while, at infinite tempera¬ 
ture, the notion of a gap still holds—entropic considera¬ 
tions dominate and the relevant gap is approximately the 
gap associated with the largest few total-spin states. The 
form implicitly assumes a dynamical exponent z = 1//?; 
the scaling of the spin size with the exponent —13/2 oc¬ 
curs because it scales as the square-root of the cluster 
size n (as long as the chain is not purely F or AF); n 
scales with the exponent —/3) at both high and low tem¬ 
peratures. The reason that both F and AF spins scale 
in the same way is attributed to the fact that proportion 
of F and AF bonds in the limit that the RG converges is 
a finite, intermediate value. Additional requirements of 
finiteness of (A) /Aq and normalization fix the form of 
the distribution. 

The emergence of a finite dynamical exponent is ex¬ 
plained as follows: in k RG levels, ^ 2^ microscopic 
spins are combined into a single cluster, while the en¬ 
ergy gap, expected to reduce by a factor r at each level, 
reduces from Aq to Ao/r^; the dynamical exponent z is 
then readily found to be z = In r/In 2. Note that even 
though our simulations with more singular initial bond- 
distributions exhibit non-universal values of the exponent 
z = 1/(3, the convergence of the RG towards a fixed value 
of the final bias indicates that the notion of a scaling dis¬ 
tribution as discussed above still applies. 

B. Calculation of Sq{uj) 

From the scaling forms of the distribution functions 
Pf(af), the dynamical structure factor Sq{uj) is given by 

SqiLu)= f dAodSLdSR6iio-Ao)N{Ao) 

a^F.AF^ 

X M^Sl, SR,)F{qno)PaiA = Aq, Sl, Sr) (2) 

The above integral simply reflects the discussion of the 
noise calculation procedure: at every step of the RG, a 
local integral of motion combing spins Sl and Sr, with 


a gap A = Aq is eliminated, producing noise of magni¬ 
tude M{Sl, SR)F{qno) at frequency oj = Aq. Here, as 
mentioned before, M{Sl,Sr) corresponds to matrix ele¬ 
ments of the spin operators associated with the transition 
at frequency w = Aq, while F{qno) modifies this ampli¬ 
tude depending on how the probe couples to the cluster 
generating the noise; this depends on the ratio of the 
size no(Ao) of the typical clusters at maximum gap Ag 
and the probe wave-length 2TT/q. N{Ao) ^ l/ng ^ Aq 
denotes the number of bonds that remain at the cut-off 
scale Aq. 

We now use Eq. (2) to evaluate the form of the struc¬ 
ture factor and show that it reproduces the results in 
Eq. (1). The noise magnitude M{Sl,Sr) ~ and, 

consequently, at any finite temperature (as any generic 
combination of spins Sl and Sr) scales as the square 
of the typical spin size, i.e., M ~ Sg ~ However, 

at zero temperature, the scaling changes to M ~ Sg ~ 
Aq ' (see Appendix A). 

Next we consider the scaling of the form factor F{qno). 
To arrive at the results in Eq. (1) starting from Eq. (2), 
we need to show that at both T = 0 and T = oo, and 
frequencies w 9 ^^^, F{qno) ~ 9 ^Uq, while for fre¬ 
quencies uj 9 ^^^, at r = 0, F{qno) ^ const, and at 
T = 00 , F^q^uo) ^ l/ 9 nQ. Let us note that the distinc¬ 
tion between the two regimes comes from the fact that at 
high-frequencies, the clusters are smaller than the typi¬ 
cal length 1/9 of the modulations of the probe, while at 
lower-frequencies, the clusters are larger in comparison. 
So, the noise generated by two clusters A and B much 
smaller than the probe scale I /9 comes with an form- 
factor F{qno) that scales as the square of the gradient 
9 of the probe, that is, F{qno) ~ q'^riQ- (Recall that 
the form factor is the square of the difference of the cou¬ 
plings of these individual clusters A and B to the external 
probe.) 

We now consider the scaling of F{qno) at low frequen¬ 
cies {uj ^ q^^d) where clusters are much larger than the 
probe scale I/ 9 . At low temperatures, clusters are com¬ 
posed of a large number of singlets, and consequently, 
cluster sizes and the bond lengths scale in the same way 
(as ^ no, see Ref. [40]). When two such clusters are 
merged, the relative phases of the probe coupling are 
effectively scrambled. Thus, in this regime, at low tem¬ 
peratures, F(qno) is 9 -independent. 

AtT = 00 , we find the scaling of F(qno) ~ l/ 9 n.g using 
the following two facts—i) singlets are entropically un¬ 
likely at high temperatures; clusters are compact (nearly 
contiguous arrays of spins that are not part of singlets) 
and bond-lengths are of the order of the microscopic scale 
and; ii) spin-clusters point in arbitrary directions; as a 
result, one can show that variance cr(g+) of the probe 
coupling of a cluster composed of two smaller clusters 
with couplings gA and gs respectively is precisely equal 
to the probe form factor associated with combining these 
clusters (see Appendix B). The first of these implies that 
the mean of the probe coupling gA of cluster A of size 
rig can be approximated by |( 9 a)| ~ /g dx cos{qx)/n ~ 
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^jqno for qriQ ^ 1. The second result shows that the 
scaling of the probe form factor is given by the scal¬ 
ing of the variance of the probe couplings; F(gno) 

= {{9A - gB f/^) ~ [ct( 5 a) + cr( 5 B)]/ 4 , where, to 
obtain the last (approximate) result, we assumed that 
the squared-mean {gA) ~ l/(gno)^ is negligible in com¬ 
parison to the variance <T{gA) in the limit qn^ 3> 1. If we 
reasonably assume that clusters A and B of similar size 
^ uq (and consequently, similar variance) are combined 
typically, then the this immediately leads to the result 
F{qno) ~ cr(g+) ~ l/^u-o (which also justifies neglecting 
the mean values {g a) jigs)) as mentioned above. 

With the aid of the specific form of the distribution 
function, the noise amplitude and the form factors dis¬ 
cussed above, we can perform the integration in Eq. (2) 
by first homogenizing all flow parameters in favor of the 
maximum gap Aq, which can then be integrated directly 
to yield the q,uj dependent behavior. This yields the re¬ 
sults in Eq. (1). 


V. GENERALIZED DIFFUSION AND THE 
STRUCTURE FACTOR 


The derivation of the structure factor presented in the 
previous section relies on specific details of the model. 
We now present a general ansatz for the spin-diffusion 
propagator in a system exhibiting anomalous diffusion 
and show that it reproduces the form of the structure fac¬ 
tor in Eq. (1) at high temperatures. (See also Ref. [15].) 


A. Ansatz for the anomalous diffusion propagator 

The spin propagator describes the decay of the spin 
density S{x,t) at any point x and time t > t' given the 
spin density profile at all points x' at some fixed t', that 
is, S{x, t) = J dx'G(x — x',t — t')S{x', t'), where, by def¬ 
inition, G{x — x',Q) = 5{x — x'). Note that we do not 
keep the tensor structure of the Green’s function because 
off-diagonal components vanish by symmetry and there 
is no preferred ordering direction in one-dimension for 
any non-zero concentration of AF bonds at any temper¬ 
ature. (In higher dimensions, the approach is directly 
applicable above any ordering temperature.) In the case 
of regular diffusion, the Green’s function satisfies the dif¬ 
fusion equation, and, in the Laplace domain, is given by 
Gq{uj) = + Dq^). This diffusive propagator has 

two salient features: i) lim(^_,,o(—iw)Gq=o(w) = 1 imply¬ 
ing that the q = 0 mode does not relax (that is, total 
spin/density is conserved), and ii) it has a pole at finite 
imaginary frequencies for all finite-q modes, which as a 
consequence, relax exponentially in time. 

We now generalize this propagator to the case of 
anomalous diffusion, where the system exhibits an 
anomalous scaling (the dynamical exponent z = 1//3 ^ 


1 , 2 ) between q,(jj: q - 

Gq{uj) = \l\-ibj F g^/^/(a;/g^/^)], (3) 

where / is some well-behaved function which ensures con¬ 
servation of the total spin [Gq^o)^) 


B. Calculation of the structure factor at T = oo 


We would like to use the Green’s function postulated 
above to calculate the structure factor Sq{ijj). First, 
we introduce the dynamical susceptibility Xg(w) which 
is the usual Kubo-response of the spin-density Sq = 
J dxS{x)e^‘>^ to a field hq that couples to S-q. That 
is, Xij('^) is the Laplace-transform of Xg(t —t') = — 


t'){ Sq{t).,S.q{t') 


The dynamical susceptibility 


is connected to the propagator Gq{uj) via the relation 
Xq{^) = X°[l + i^Gqiuj)], where x° is the static suscep¬ 
tibility at wave-vector q. The result can be rationalized 
as follows: the measurement of the dynamical suscepti¬ 
bility is carried out by slowly ramping up the perturba¬ 
tion hq{t) = (coupled to S-q) for times t < 0 and 

observing the relaxation of the spin-density Sq for subse¬ 
quent times. One can think of such an experiment as one 
that sets up a spin-density Xq^q by time t = 0, and whose 
subsequent relaxation is given by the diffusion propaga¬ 
tor, i.e., {Sq{u!)) = XqhqGq{u!) (ill Laplace-domain; note 
that {Sq{u!)) denotes the expectation value of the spin- 
density operator and not the structure factor). Alterna¬ 
tively, we can appeal to the definition of the dynamical 
susceptibility to directly find {Sq{u})) = [x9(w)—Xq]Iig/*w 
(see Sec. (7.14) of Ref. [43]). These alternative interpre¬ 
tations yield the relation between Gq{uj) and Xq{'-^)- Fi¬ 
nally, the structure factor is related to the imaginary part 
of the dynamical susceptibility using the (fluctuation- 
dissipation) relation Sq{u)) = coth (tLi/2T)Im [xq{uj)]. 

At high temperatures (w <C T), the Green’s func¬ 
tion in Eq. (3) yields a structure factor Sq{u!) = 
Tx\q-^l^f{^/q^lf^)/\{.^lq^l^? + It is easy 

to confirm that the above form can be represented as 
Sq{LS) = const. X g“^/^5'oo(w/g^/^), in agreement with 
the result for S^{oj) in Eq. (1) if we assume that Tx° 
has a finite limit for small q <C 1. This is expected for 
any physical system with short range interactions, at high 
temperatures. 

To complete the argument, we now determine the lim¬ 
iting forms of the scaling function gaoiv = 
that is, its high (low)-frequency limits y 1 {y 1). 
First, note that the (real part of the) optical conduc¬ 
tivity CFq{u}) = {[Jqiuj)., J-q(—u})]) /lu) cuu be computed 
from the structure factor using the relation Re [cfq{ui)] = 

tanh (w/2T)5'q(a;); the result (at any non-zero fre¬ 
quency) follows from the continuity relation iuiJq{uj) — 
iqSq{ui) = 0 and using the fluctuation-dissipation re¬ 
lation discussed above. Straightforward calculation us¬ 
ing these results yields a{q = 0,w) = c{q = 0)a;^“^^, 
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where the pre-factor c{q) is finite and non-zero only if 
9oo{y) ~ for y ^ 1 . Since there is no reason 

for the optical conductivity to vanish or diverge at fi¬ 
nite frequencies in a high-temperature system, we expect 
this condition yields the scaling of the function gao{y) 
in the high-frequency limit. The argument for the low- 
frequency limit is as follows. The form of the Green’s 
function ensures G{x = 0,t) ~ 1/t^. On this basis, 
we expect large q modes should decay algebraically as 
~ 1/t^ as well. This implies that /(y ^ 1) ~ y^“^, 
and consequently, gaoiv ^ 1) . These results 

together reproduce the complete form and limits of the 
structure factor (T = oo) in Eq. (1). Because of 

the generality of this derivation, we expect that any sys¬ 
tem exhibiting anomalous diffusion must have a structure 
factor that scales as described in Eq. (1). 


C. Failure of linear response at T = 0 

The generalized-diffusion ansatz fails to explain the 
T — 0 structure factor. In particular, the low-frequency 
power-law Sq{u!) ^ l/w“, with 0 < a < 1 is impos¬ 
sible to obtain with this Green’s function form. Fur¬ 
thermore, according to the thermodynamic sum rule, 
Xiq) = /dwx"(w)/w (x"(y,a;) = sign(w)S',(a;)), the 
static susceptibility is expected to diverge even at finite 
q making the derivation suspect. Hence, at T = 0, we 
are forced to conclude that either linear response does 
not work, or the scaling hypothesis for the Green’s func¬ 
tion breaks down, or both. In what follows, we provide a 
heuristic argument for the breakdown of linear response 
by showing the the divergence of the static susceptibility 
at finite-y. 

Following Ref. [40], we imagine applying a perturb¬ 
ing field to the chain. This perturbation cuts 

off the RG flow at a scale where the energy gap Aq is 
comparable to the Zeeman energy of typical, polarizable 
(of size n < 1/q) spin-clusters; that is, the RG is cut¬ 
off at Aq ~ hq min [i/ho, 1 /y^] , where we note that 
the typical cluster size at an energy scale Ag is given 
by no ^ 0 ^: and possesses a spin ~ i/ho- Grucially, 
the magnitude hq determines the point at which the RG 
is cut-off (in particular, whether no > 1 /y or ng < 1/q 
at this scale) and consequently, the response of the sys¬ 
tem. Let us define no{hq) and Ag(/iq) as the amplitude 
/iq-dependent length and energy scales at which the RG 
is cut-off. The condition no{hq) ^ 1/q can be alterna¬ 
tively cast (using the scaling results) as q^G+'^/P ^ 
and both limits of this condition may be experimentally 
relevant. 

We first analyze the case q ^ l/no{hq), Ao{hq) ~ 
hq/ y/q first. We note that the magnetization ni of polar¬ 
ized clusters of size n is yjn for n < 1/q and small oth¬ 
erwise. If we assume that the distribution D{n) of the 
cluster size n satisfies D{n < ng) ^ I/ug and D{n > no) 
is negligible (the form of D(n) is numerically justified, see 
Appendix G), then the average magnetization is given by 


ruq = nQ^(hq)q^'^^ = q hh/P. Thus, for 2/? < I, the 
static susceptibility Xq = nT-q/hqlh^o diverges even at 
finite q. 

For the opposite limit of y ^ l/no{hq), Ao{hq) ^ 
hqy/n/, the magnetization is simply niq = m^^g ~ 
l/\fno and the susceptibility diverges for all values as 
P as, Xq = mq/hq\h^o = in agreement with 

the y = 0 result of Ref. [40] . 

Thus, the static susceptibility generically diverges at 
T = 0 in this system, although the precise nature of the 
divergence is controlled by the condition q^G+^/P ^ 

Let us note that, a divergent static susceptibility at cer¬ 
tain wave-vectors for clean F (at y = 0) or AF (at 
y = tt) systems is not surprising—a small perturba¬ 
tion with the correct wave-vector results in a macro¬ 
scopic reduction in free-energy and consequently, a di¬ 
vergent response. Since, in our disorder system, the 
probability distribution of cluster sizes, D{n), has sig¬ 
nificant weight for all clusters of size n < uq (and ng 
diverges in the zero-frequency limit), our system has a 
thermodynamically-relevant presence of approximately 
independent spin clusters of all sizes. This results in a 
divergent response to perturbations over a broad range 
of wave-vectors. 


D. Relation of 1/f noise and 
many-body-localization at finite temperatures 

For frequencies uj /2^“ (where zg is the dynamical 
exponent found for the system at T = 0, and is non-finite 
only in the purely AF case), we expect the behavior of the 
structure factor to be given by the infinite-temperature 
form of the structure factor we compute. This structure 
factor has a noise spectrum l//“, where a < 1. We can 
now reverse the arguments of the previous section to con¬ 
clude that the spin-propagator G{q,u}) ^ l/w°‘ for fre¬ 
quencies u) <C T/2^° ,q^/P where j3 = 1 — a. Fourier trans¬ 
forming this result implies directly that G(y, t) ~ 
that is, the spin-density at wave-vector y decays alge¬ 
braically with a power 1 — a > 0 , where, in particular, if 
a = 1, the mode does not decay completely even at in¬ 
finite time. The inability of certain non-conserved quan¬ 
tities to decay to equilibrium values is a defining char¬ 
acteristic of many-body-localization, and it is reasonable 
to posit, on the basis of this discussion, that when the 
noise power law is precisely a = 1 (or, noise is 1 /w), the 
system is at the verge of becoming many-body-localized. 

A complimentary picture is obtained by looking at the 
dynamical exponent I//3 = 1/(1 —a) (at finite T). In the 
course of RG, when a cluster joins another cluster to form 
a supercluster, its spin begins to process around the di¬ 
rection of the total spin. Further mergers into ever larger 
clusters add new axes of precession with ever decreasing 
frequencies, which leads to scrambling of the spin polar¬ 
ization. The RG suggests that, for the disordered Heisen¬ 
berg system, the time scale for the phase scrambling over 
a length scale I scales as a power law ti ~ l^^P. This im- 
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plies that for a < 1, it takes only polynomial time in the 
distance to transport spin across this distance, while, for 
a = 1, it becomes exponential (or worse) in the distance, 
which is suggestive of MBL. 

Since we do not observe such a case (with a = 1) in our 
numerical simulations, we conclude that the disordered 
Heisenberg spin chain is not localized. Nevertheless, we 
note that as the disorder strength is increased, the noise 
exponent a tends towards 1. Moreover, at low temper¬ 
atures, the power law a = 1 is obtained for the purely 
AF chain which flows towards the m/in*te-randomness 
point—here the integrals of motion are exact since the 
RSRG is precise in the low-energy limit. 


VI. SIMULATIONS IN TWO-DIMENSIONS 

Simulations in two-dimensions are performed analo¬ 
gously to the one-dimensional case except that the num¬ 
ber of neighbors of an effective spins can change over the 
course of the RG. Maintaing this extended network in¬ 
troduces a computational cost that restricts simulations 
to smaller system-sizes in two dimensions; we consider 
square systems of spins on a 50 x 50 lattice and ‘strips’ 
of spins on a 500 x 6 lattice. The growth in the connec¬ 
tivity of the network over RG steps also reduces the effi¬ 
cacy of the RG because achieving a separation of scales 
becomes less probable. In particular, at T = 0, when 
singlet-formation is preferred, the clusters become highly 
non-compact, extensive objects with a large number of 
neighbors, and consequently the fidelity, which is the ra¬ 
tio of maximum energy scale in the system compared 
with the energy scales amongst the spins neighboring the 
maximally-coupled spins, does not rise substantially over 
the course of the RG procedure, even if the distribution of 
couplings assumes a power law form (see Appendix I). At 
high-temperatures, the RG is more robust since singlets 
are entropically unfavorable, and clusters grow uniformly 
keeping the number of neighbors relatively unchanged 
over the course of the RG. In what follows, we restrict 
ourselves to a discussion of the T = oo simulations on 
the two-dimensional square and strip geometries. 

For a spin network arranged in a square geometry, we 
find that the structure factor exhibits a scaling collapse 
in accordance with the generalized-diffusion form [S'^{uj) 
in Eq. (1)] with a dynamical exponent z « 2; this is ac¬ 
companied by a low-frequency tail in the finite-g 

structure factor (see Appendix I). We note that the dy¬ 
namical exponent z = 2 is also associated with regular 
diffusion; the important difference here is that the diffu¬ 
sive structure factor is flat at low frequencies. Since the 
fidelity does not improve significantly over the course of 
the RG for systems with square geometries, we cannot 
conclude definitively on the nature of dynamics of such 
systems. 

From the point-of-view of experiments measuring flux 
noise in SQUID systems, we consider the case of a two- 
dimensional strip of spins (the width of the conducting 


strip of the SQUID circuit is typically much smaller than 
the length of the SQUID loop). In this case, the RSRG 
finds that the structure factor possesses a 1/v^ power- 
law frequency-dependent behavior which eventually, at 
lower frequencies, gives way to an anomalous power-law 
behavior, l/w“; a « 0.70 (see Fig. 16). The reason 
for this is self-evident; the clusters contributing below 
a certain frequency are wider than the width of the spin 
network and can be thought to be lined up in an effec¬ 
tively one-dimensional array; thus, the inheritance of the 
anomalous power-law behavior from one-dimension is un¬ 
surprising. 


VII. SUMMARY AND OUTLOOK 

In this work, we extended the RSRG protocol with 
rules for the renormalization of the couplings of a generic 
probe of spin fluctuations. This allowed us to compute 
noise spectrum measured by such a probe, and in partic¬ 
ular, the dynamical structure factor Sq{uj). We find, for 
one- and quasi-one-dimensional systems with exchange¬ 
coupling magnitudes distributed according to a 1/| J| dis¬ 
tribution, and a wide range of initial composition (pro¬ 
portion) of F/AF bonds, that the RG converges at both 
T = oo and T = 0 to non-universal strong-randomness 
fixed points associated with an anomalous dynamical ex¬ 
ponent z. This dynamical exponent is associated with 
the spectral characteristics of the noise power a at low 
frequencies: a = l — 1 /zatT = oo and a = 1 — l/2z 
at T = 0. We also find a non-trivial scaling of the high 
frequency end of the structure factor [see Eq. (1)] and 
relate the power law of the frequency dependence in this 
regime to the low frequency noise power a. 

These findings were justified for both high and low 
temperatures by postulating a scaling form of the prob¬ 
ability distribution function Pf{af) governing the distri¬ 
bution of F(AF) bonds with a given coupling, and at en¬ 
ergy scales where the RG has converged. Additionally, we 
showed how an ansatz that generalizes the spin-diffusion 
propagator, reproduces the complete form of the struc¬ 
ture factor at high temperatures, from which we conclude 
that the results obtained (including the presence of l/a;“ 
noise) are universal and extend generally to systems ex¬ 
hibiting anomalous diffusion. We explain that the failure 
of this approach in the zero-temperature case occurs due 
to the failure of linear-response theory for the random- 
bond Heisenberg spin chain at this temperature: even 
at finite wave-vectors, static susceptibility of the chain 
diverges. 

From the point of view of interest in many-body- 
localization physics, our work raises an important 
question—is the phase we find non-localized but also non- 
ergodic? (Such a phase was proposed in Ref. [39].) The 
RSRG protocol creates a macroscopic number integrals 
of motion (albeit increasingly non-local integrals) which 
implies the breakdown of ergodicity. However, since 
the RSRG protocol suggests a flow towards a strong- 
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randomness fixed point, and not m/inite-randomness, the 
RG protocol is not exact and neither are the integrals of 
motion it creates. This does not exclude the possibil¬ 
ity that integrals of motion do exist, and one needs to 
devise a more elaborate form of the RSRG protocol to 
correctly capture them. Thus, the status of ergodicity in 
our model is unclear. 

Irrespective of the resolution of this discussion, let us 
note that the system is certainly not in a many-body- 
localized phase. As was mentioned above, the l/a;“ 
{a < 1) form of the finite-g structure factor implies that 
all spin-waves decay in time as As long as a < 1 

(strictly), we expect polarization in the system to relax 
(to equilibrium). A complimentary picture is obtained 
by looking at the dynamical exponent 1//3 = 1/(1 — a) 
(at finite T) which governs the scaling between the time 
t it takes for spin phase to get scrambled over a distance 
1; t . As a —)■ 1, the scaling changes from polyno¬ 

mial to exponential (or worse) in distance, which again 
suggests that the localization transition occurs at a = 1. 
Since we never observe this value for a in the infinite- 
temperature simulations (although, increasing disorder 
strength D is seen to slowly drive the exponent a to¬ 
wards 1), we conclude that the disordered Heisenberg 
spin chain is not localized at finite temperatures. 

Finally, we note that the noise generated by these 
random-bond Heisenberg spin chains with a 1/|T| distri¬ 
bution captures some interesting facets (as seen in Fig. 4 
for the experimentally relevant case of equally propor¬ 
tioned F/AF system, see also Appendix D) of exper¬ 
iments observing flux noise in SQUIDs: i) the noise- 
exponent decreases with increasing temperature; ii) the 
magnitude of the noise is fairly temperature indepen¬ 
dent over a large frequency range; and hi) the point 
of constancy in the magnitude of noise spectrum (as a 
function of temperature) is about 9 orders of magnitude 
smaller than the microscopic energy scales. The caveat 
is, of course, that the network of spins in SQUIDs is 
two-dimensional and these simulations were performed 
for spin chains—due to computational costs, we can¬ 
not perform simulations on equally large two-dimensional 


systems and are, consequently, limited in the frequency 
range of the calculated noise. Nevertheless, the simu¬ 
lations performed on elongated two-dimensional strips 
of spins show, below a certain frequency, an anomalous 
power-law frequency-dependence akin to that seen in one¬ 
dimensional chains. We believe that this may be a rea¬ 
sonable representation of the experimental setup. 

In future work, we aim to extend our RSRG protocol to 
systems with anisotropic dipolar couplings and (or) ran¬ 
dom local magnetic fields that create an additional pre¬ 
ferred axis amongst the spins; such couplings are known 
to change the dynamics of system significantly [44, 45], 
making it resemble the dynamics of disordered Ising 
spins, which furthermore, are likely to be more amenable 
to a RSRG analysis even in higher dimensions. Moreover, 
in one-dimension, such models are known to exhibit the 
many-body localization transition [10, 46] and will thus 
allow us to further probe the connection of l/w noise with 
the many-body-localization transition. It would also be 
interesting to examine the presence or absence of l/a;“ 
noise in other systems where anomalous dynamical ex¬ 
ponents are observed; such as non-Fermi liquids [47, 48] 
and Luttinger liquids [49]. 
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Appendix A: Calculation of Noise 


To the zeroth-order in the perturbative parameter 
A/Ao (A is the gap between successive energy levels of 
the spin-pairs formed due to the neighbors of the spins in 
the pair forming the maximum gap Aq), the noise con¬ 
tribution at every step of the RG is given by: 




+ 9 


A 


\ ^ 

s — l,rn') 6(uj J- w^) 
s-|-1,to'^ (5(a; —(A1) 


where A = S- — 


S-.S+ 


|S+| 


^>5'+, 9^ = {9A - 5 b )^/4 is the 


geometrical factor proportional to the square of the dif¬ 
ference of the coupling of the probe to the two spins A 
and B being combined at the given RG step, and a;^( 
is the energy difference between states of total angu¬ 
lar momentum s and s J- 1 (s — 1) formed from the two 
spins. We will equivalently refer to these as the frequen¬ 
cies of the ‘up’ and ‘down’ transitions associated with 
a state with total momentum s. m and m' are the az¬ 
imuthal quantum numbers associated with the states s 
and s J- 1 (or s — 1 in the second line). Note that there 
is no sum for to: all —s < m < s contribute equivalently. 
Moreover, at this step of the RG, we only choose the 
magnitude s of the effective spin, and not its completely 
quantum state. 

The spectrum of two spins is readily given in terms 
of the coupling strength J between the spins of size sa 
and sb being combined; the energy for such a spin pair 
combining into an effective spin of size s is given by 
Es = i (s(s + 1) - sa(sa + 1) - sb(sb + !))• The en¬ 
ergy difference between states s and s ± 1 can be com¬ 
puted readily from this result. It is important to note 
that the sign of the Ferromagnetic and AntiFerromag- 
netic coupling implies a certain ordering of these levels, 
which changes the meaning of and accordingly. 

S- = Sa — Sb is a rank 1 spherical tensor operator (as 
one can ascertain from its commutation relations with 
the total spin angular momentum operator 5'+). There¬ 
fore, it only leads to transitions between states s and 
s ± 1 , besides also having a matrix element that leaves 
s unchanged. The projection of S- on to S+ (which is 
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subtracted from S- to yield A) precisely removes this 
matrix element corresponding to transitions from s —>■ s; 
this projection operator is a low-frequency component of 
the two spins that must not be eliminated at the given 
RG step. 

The remaining aim of this sub-section 

is the calculation of the matrix element 


E. 


s ± 1, m' > ( s ± 1 , to' 


A 


By 


Wigner-Eckart theorem, tbe matrix elements are related 
to the Clebsch-Gordan coefficients as 


(s, m\ Aq |s', to') = Cs,s' < s, to||s', to'; 1, g > . 

Thus, the up-transition rate is given by 

m' ,q 

and similar expression for M{s,].). Wigner-Eckart does 
not provide the values of the coefficients Cs^s >; they can 
be established, however in the following way. Gon- 
sider the matrix elements of Aq, (s, to| Aq |s', to') = 
Cs,s' < s, to||s', to'; 1 ,0 > and (s', to'| Aq |s, to) = Cs\s < 
s', to'I |s, to; 1,0 > . They are obviously complex conju¬ 
gates of each other. Thus, 


_ ^ < s', to'||s, to; 1,0 > _ * 2s'-I-1 
^ < s,ui||?;TO'rb 0 > ^ 2EEr’ 

where the last equality follows from the properties of 
symbols. On the other hand, 

M (s) = M (s, j') -f iiT(s, j,) =2si(si-|-1)-|-2s2(s2-|-1) 

( , (si(si + 1) - S2(S2 -I- 1))^ 

- s(s 4-1)-p-—-. 

^ ’ s(s-f 1) 

(A2) 

These two relationships are sufficient to calculate all tran¬ 
sition rates separately. Eor the purposes of Sec. IV B, it 
suffices to note that, generically, the sum M(s) scales as 
~ s^, however, for s = si -I- S2 and s = si — S2, i.e., 
extremal states, as is the case in the T = 0 RG flow, 
M(s = Si -I- S2) = 4 siS 2 /(si -I- S 2) and M(s = si — S2) = 
4(si -I- l)s2/(si - S2 -f 1), implying, at T = 0, M(s) s. 


Appendix B: Analytical argument for scaling of the 
form factor at T = 00 

In the main text, we argued that the probe form factor 
E(gno) ~ 1/guo for qriQ ^ I in the infinite temperature 
case. Below we provide a more detailed justification of 
this result. 

Let us first recall that the probe form factor F{qno) 
describes the scaling form of the squared difference of 
couplings gA and gs, that is, — gs)^I'i, when these 
individual couplings correspond to clusters whose size 


is ~ riQ (note that the probe couplings contain the q- 
dependence). In this appendix, we show that this factor 
is related to the variance in the distribution of probe cou¬ 
plings which scales as l/guo for qn^ ^ 1. 

The general result for the effective coupling g+ of the 
probe to a cluster formed by merging two clusters A and 
B of spins Sa and Sb, coupled to the probe by strengths 
gA and gs, is given by 


9+ = 


gA + gs , gA- gs 




Sb 


S1 + S2 


(Bl) 


This expression in Eq. (Bl) is valid for quantum spins; 
however, at the advanced stages of RG at T = 00 we can 
treat them as classical variables. Then by averaging uni¬ 
formly over spin directions, we can evaluate the average 
i9+)sA,SB = ^9a + gs)/^ and and the variance of g+, 

<^'^{9^)sa,Sb = {gVlsA^Sn ~ <5+)L.Sb = {9A-9B?I^- 

Here (--Ea Sb ioaplies averaging over spin orientations 

of the 2 constituent spins Sa-, Sb of a cluster; note that 
gA and gB will themselves have a distribution of values 
owing to the fact that they are themselves composed of 
smaller spin clusters. In what follows, we will use the 
notation (...) (without the subscripts) to imply the aver¬ 
aging over all internal spins of a cluster. Erom the above, 
we see that the noise factor F{q£) obtained when com¬ 
bining two clusters of size i with probe couplings gA and 
gs (which depend on g) is given precisely by the variance 
cr^(g+). Thus, in order to find the noise factor F{q£), we 
need to determine the distribution of g+. 

The expression for (g+)g^ implies that the probe 
coupling to each cluster is approximately the average of 
the probe couplings to individual spins. Note that spins 
which form singlet pairs do not couple to the probe, and 
are not counted in this sum. However, as argued before, 
singlets are entropically unfavorable at T = 00 and so 
clusters are primarily arrays of contiguous spins. Gonse- 
quently, the | (g+) | ~ Jq dx cos{qx)/i. In particular, for 

g » 1/^, I (5-h) I « iM 

Now we estimate the width (tT(g_|_)) = 

((g^ -I- g^ — 2g^gB)/4). We will self-consistently 
show that the factor (g^gs) can be dropped from this 
sum if these clusters have size £ ^ 1/g. To this end, 
we first note that the phase difference between gA and 
gB is of the order of which is effectively random 
for q£ 3> 1. Thus, such clusters have uncorrelated 
probe couplings and (g^gs) « {9a) {9b) ^ l/(g0^; 
the last relation is determined from the result of (g+). 
Assuming that this term makes an unimportant con¬ 
tribution (in the large £ limit) to the variance, we find 
(tT(g+)) = ((g^) + (g^))/4. If we merge clusters with 
similar lengths £ and similar values of (g^) ~ {9b}^ 
the resulting average value (cr^(g+)) goes down by a 
factor of 2; consequently, (cr^(g+)) ^ l/£ for £ 1/g. 

In particular, since the width of distribution begins to 
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grow only for i ^ 1/q, cr‘^{g+) ^ l/?^- This scaling 
self-consistently confirms our assumption to neglect 
the cross term {gA) {gs) ~ completes the 

argument. 

Therefore, the probe form factor F{qnQ), which is the 
scaling form of the factor {gA — gB)"^ 1^ when the size of 
clusters A and S is ~ no, scales as ^ l/{qno) for quo ^ 1, 
as advertised in the main text. 





FIG. 5. (a) The T = oo dynamical structure factor (cj) 
is plotted as a function of frequency for various q = 2n/l 
{I G [40,800]), system size L = 15000. (b) S^{u)) ~ 1/q for 
uj q^^^ is evident from the scaling collapse at low frequen¬ 
cies. {c)S^{uj) ~ q^ for LJ 3> q^^^ is evident from the scaling 
collapse at high frequencies. 



10-8 iq-4 iqO 


FIG. 6. (a) The T = 0 dynamical structure factor Sq{uj) 

is plotted as a function of frequency for various q = 271-/1 
{I G [40,800]), system size L = 15000. Sg(u) ~ const, forw ^ 
is evident from the collapse at low frequencies (without 
any scaling), (b) Sq(uj) ~ q^ for lu ^ q^^^ is evident from the 
scaling collapse at high frequencies. 


Appendix C: g-dependent scaling of the Structure 
factor 

Figs. 5 and 6 detail the g—dependent behavior of the 
structure factor S//°{u}) at T = oo and Sg{uj) a.t T — 0 
respectively, at high (w » q^^^) and low (w <C q^^^) 
frequencies. The g—dependence is in accordance with 
the results in the main text. 


Appendix D: Comparison of T — 0 and T = oo noise: 


Fig. 7 shows a direct comparison of zero and infinite 
temperature noise for single disorder realization. There 
are three aspects of this plot that are of significance: i) 
the magnitude of the noise does not vary much with tem¬ 
perature, while the power laws change slightly; ii) (for the 
case of a 50-50 mix of F/AF bonds) the low frequency 
power law seems to increase with decreasing temperature, 
and; iii) the noise magnitude seems to be fairly constant 
at approximately w = 10“®. A curious observation in 
nearly all experiments is that the noise magnitude does 
not change at approximately IHz, which is about a factor 
of 10“® smaller than the typical microscopic interaction 
strength amongst the spins. In contrast, in our simula¬ 
tions, the microscopic scale is ~ 10, = 3, and noise 

magnitude is observed to be roughly constant at about 
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uj « 

We note that the above comparison is based on the 
assumption that the two-dimensional network of spins 
producing flux noise in SQUIDs is likely more tightly 
confined in one direction (along the width of the con¬ 
ducting strip) than the other (along the length of the 
conducting strip), and that this implies that the behav¬ 
ior of the spin-network mirrors that of one-dimensional 
spins below a certain frequency scale. 



FIG. 7. The structure factors S^{lo) (T = oo, blue) and 
Sg(Lu) (T = 0, red) are plotted simultaneously, for initial bias 
rji — 0, disorder strength D = 3, system size L = 15000, 
and q = 27r/160. The powers a' and a are also indicated 
for the two cases with the corresponding subscripts. They 
approximately satisfy the respective scaling relations for zero 
and infinite temperature cases. 


to the spin-pair being eliminated at the given RG step, 
and A is the maximum gap formed by the coupling of 
one of these spins (in the spin-pair) to its neighbors. It is 
seen to rise dramatically as the RG reaches convergence 
indicating that the RG is fairly successful; since the fi¬ 
delity is large but finite we conclude that the RG flows 
to a strong-randomness fixed point, and not an infinite¬ 
randomness point. In Fig. 9, the fidelity is plotted for 
the infinite temperature RG simulation and is seen to 
also perform quite well. Note that, at zero temperature, 
the fidelity is best for the initially completely AF sys¬ 
tem, and becomes progressively worse with the addition 
of more F bonds. 

Fig. 10 plots the distribution of the gap A after 98% of 
the spins were eliminated for the same system as above. 
The power law form of the distribution guarantees that 
the typical gaps are much smaller than the largest gap 
in the system, which ensures that the RG works (and 
increases the fidelity). The noise power law a and the 
power law of the distribution of the gap is not the same 
generally. The the purely AF case, the distribution is 
Da (A) ~ 1 /A, as is known from the analysis of the IRFP 
in the AF system at T = 0. As the initial distribution 
becomes more ferromagnetic, the exponent of the power 
law approaches 0. The smaller exponent also marks a 
reduction in the fidelity and the efficacy of the RG. 

The power 7 of the gap distribution Pa (A) ~ 1/A^ 
and the hnal bias 77 / are plotted as a function of the initial 
bias r]i of a 1/| J| distribution with disorder strength D = 
3 in Fig. 11. 


Appendix E: Numerical checks / convergence of RG 

Fig. 8 plots (a) bias, (b) noise frequency and (c) ‘fi¬ 
delity’ is plotted against the RG step for a system with a 
given set of initial conditions (bias, disorder strength, sys¬ 
tem size, temperature). Although plotted for a specihc 
system, the following observations are fairly independent 
of the various system parameters, including temperature. 
The bias is seen to flow to a a stable value after about 
90% of the initial spins have been eliminated. This value 
is close to the 50% - 65% (77 = 0 - 0.3) mark although it 
has small variations which depend on the initial bias of 
the system at T = 0, but not at T = 00 . The noise fre¬ 
quency (of eliminated spin-pair) dips below microscopic 
values only after these initial eliminations. These obser¬ 
vations suggest that the RG converges and the (high- and 
low-) frequency dependent noise observed at both ends of 
the crossover (w ^ is from the limit where the RG 

has converged. The fidelity plotted is the log-averaged 
(typical) estimate of the inverse of the perturbation pa¬ 
rameter, that is, / = Ao/A, where Aq is the gap due 


Appendix F: Relation between power laws of the 
gap distribution and the noise exponent 

The power laws of the gap distribution and the noise 
exponent are related. This can be seen as follows. The 
number of bonds at a given maximum gap value Aq 
is given via the dynamical exponent z as Nt(Ao) ^ Aq' . 
Some of these bonds are F and some AF; we assume that 
the proportion F/AF bonds at the maximum gap Ag is 
PF/Api^o)- We assume a power-law form of the distri¬ 
bution of the F/AF bonds associated with gaps A given 

by (3j’/af(A, Ag) = In the regime 

that the RG has converged, the exponents 'Jf/af and 
proportions Pf/af become independent of the maximum 
gap Ag. 

We now prove the relation 1/z = pf{1 — Jf) +PAFi^ — 
'Jaf) as follows. The number of F/AF bonds elimi¬ 
nated at every step of the RG is given by dNp/AFi^o) = 
At(Ag)(5p'/^p'(Ag, Ag)dAg. The sum of the eliminated 
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FIG. 8 . The average (a) bias, (b) noise frequency and (c) 
‘fidelity’ is plotted against the RG step for a system of size 
L = 10000. Initial bias r]i = 0.5, temperature T — 0, disorder 
strength D = 1. 



FIG. 9. The fidelity is plotted as a function of RG steps at 
T = oo for a system of size L = 15000 and initial disorder 
strength D = 1. 


F/AF bonds gives the number of total eliminated bonds, 
which can alternatively be determined through the dy¬ 
namical exponent as (iA^t(Ao) = ^ compari¬ 

son of the results obtained yields the aforementioned 
relation. Note also that the fact that the propor¬ 
tion of F/AF bonds becomes constant implies that 
dNp{Ao)/dNAF{^o) = Pf/paf- This further implies 
that 'yp = jaf- 

The above results can also be used to determined the 



FIG. 10. Distribution of the gap A after 98% of the spins were 
eliminated. The system has Initial bias rji = 0.5, temperature 
T = 0, disorder strength D = 1. The noise power law a and 
the power law of the distribution of the gap are not the same. 



FIG. 11. The power 7 (red) of the gap distribution Pa(A) ~ 
l/A'*' and the final bias 7 / (blue) are plotted as a function of 
the initial bias of a 1/| J| distribution with disorder strength 
D = 3. 

noise exponent via the usual relations a = 1 — l/2z (at 
T = 0) and a = 1 — 1/z (at T = 00). These relations 
are discussed in Sec. IV In particular, for the case of the 
disordered Heisenberg AF chain, we know that pp = 0, 
PAP = 1, JAF = 1, and consequently, a = 1. 


Appendix G: Distribution D{n) of cluster size at 

T = 0 


It was mentioned in the main text that the distribution 
of clusters at T = 0 follows a form D{n) ^ l/uo for 
n < uq and small for n ^ uq, where ng is the typical 
cluster size at any RG scale. Such a distribution was 
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FIG. 12. The distribution D{n) of the cluster sizes n at var¬ 
ious (different colors) stages of the RG for L = 15000 spin 
chain at T = 0, and initial bias r]i = 0.5. The features of 
the result are not dependent on the specifics of the system 
parameters. The important facet here is that D(n) has ap¬ 
preciable weight at all stages of the RG in the n —>■ 0, while it 
decays beyond a certain scale no which depends on the energy 
scale at which the RG was terminated. This justifies the form 
of the distribution we use in the main text to calculate the 
T = 0 finite-Q susceptibility of the spin chain. 

shown to lead to the failure of linear response theory 
at T = 0, primarily because it guarantees the presence 
of clusters of small sizes (and especially, of the probe 
wavelength lulq) at all energy scales of the RG; these 
small clusters were found to give a divergent contribution 
to the susceptibility of the system even at finite-g. In 
Fig. 12 we justify this form of the distribution function 
numerically. 


Appendix H: Correlations between the coupling and 
the spin size 


It was argued in Ref. [27] that the RG procedure fails 
because it generates exceedingly large spins. This conclu¬ 
sion is incorrect; it is important to note that these larger 
spins also couple much more weakly to their neighbors, 
as seen in Fig. 13, and that this balance precisely en¬ 
sures that the RG can flow to a strong-randomess fixed 
point with tractable scaling properties. The authors of 
Ref. [27] perform RSRG by keeping the length of the 
chain fixed and attaching new spins with arbitrary sizes 
but also arbitrary coupling strengths which misses out 
on these correlations. Moreover, they perform their RG 
by eliminating the largest bond value J rather than the 
larger gap A. 


Appendix I: Simulations in two dimensions 



FIG. 13. Scatter-plot of the logarithm of the couplings log{J) 
of effective spins of size s remaining after the RG procedure 
has reached convergence operating on a L = 5000 spin chain 
at T = 00 . The maximum coupling is seen to decrease with 
the spin size s 


The basic features of the results in 2D (at T = oo) are 
outlined in Fig. 14. The simulations (square geometry) 
were carried out for a 50 x 50 lattice of spins. The 
fidelity does not improve markedly over the course of the 
RG. However, the computed structure factor shows ap¬ 
proximate scaling collapse according to the Generalized- 
diffusion form of the structure factor with an inverse dy¬ 
namical exponent /3 « 0.55 ± 0.05. 

While the RG doesn’t successfully converge for large 
systems with a square geometry, it still manages to re¬ 
cover the full energy spectrum for sufficiently disordered, 
small systems with reasonable success (see Fig. 15). To 
simulate a two-dimension strip of spins, simulations were 
carried out for a 500 x 6 lattice of spins. Below a cer¬ 
tain frequency scale, clusters become of the size of the 
width of the network, and the l/a;“ noise behavior of 
one-dimensional systems with an anomalous exponent a 
is recovered (see Fig. 16). For a disorder strength of 
D = 3, a Ki 0.7 was found. 

Appendix J: Role of faraway resonances 

Here we present an argument (based on a generaliza¬ 
tion of the analysis in the supplementary of Ref. [50]) 
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FIG. 14. Various features of the simulations on two- 
dimensional system of size 50 x 50 are shown. (a) 
The noise spectrum for low frequencies has a power law 
tail Sq{uj) ~ {q = here); (b) The distribu¬ 

tion of gaps is fairly broad; (c) The fidelity / does not 
rise over the course of RG; (d) Approximate scaling col¬ 
lapse using the Generalized-diffusion ansatz using exponent 
P = 0.55±0.05 for wave-vectors q — ( 27 rn/L, 27mlL) with n G 
[2(blue), 4(red), 6(magenta), 8(cyan), lO(black), 12(green)]. 
and L = 50. The lower (upper) frequency power laws 
are found to be: a ~ 0.5 ± 0.05, cx « 1.75 ± 0.05; it is 
unclear whether there is data for a sufficiently large range of 
frequencies to extract the high-frequency power law reliably, 
especially since it exhibits a drift to higher values as we go 
to smaller wave-vectors. 


that suggests that once the system develops an dynam¬ 
ical exponent z = 1//3 > 1, faraway resonances can 
be neglected. On the basis of the numerical evidence, 
we assume a power-law distribution of gaps D{A) ^ 
(^) ■ examine the possibility that two far¬ 

away spin-pairs, with local gaps separated by dA be¬ 
come resonant with one-another. The density of spin- 
pairs separated by a local gap dA is p = D{AQ)dA — 
(1 — 7 )(c?A/Ao) and the typical distance between such 
approximately-resonant spin-pairs is I = 1/p. We would 
like to compare dAo to the effective coupling between 
these spin-pairs. Since there is a fixed dynamical expo¬ 
nent 1//3 in our system, this effective interaction between 
the spin-pairs is given by Aint = , where I is again 

the typical distance found above. 

For faraway resonances to be unimportant, we require 



FIG. 15. The energy spectrum of a small 3x3 lattice spin 
system (disorder strength D = 3) is calculated using exact 
diagonalization (blue) and compared with the RG generated 
spectrum (red). 


dA ^ Aint{dA) = const. (dA)^/^. This condition is 
always satisfied if /3 < 1, a condition which is fulfilled 
in all our simulations. Thus, we expect resonances to 
be unimportant, so long as the RG generates a power- 
law distribution of gaps, and the dynamical exponent 

z = 1/^ > 1. 


Appendix K: Holstein-PrimakofF method for T — 0 

The RG seems to generally fail for a system that is pri¬ 
marily ferromagnetic. This is indicated by consistently 
poor fidelity (/ ~ 2) and a distribution of gaps that fails 
to develop a singular form. Thus, we need an alterna¬ 
tive method to understand the dynamics in the nearly 
F spin chain. To do so, we analyze the purely F chain 
by performing a Holstein-PrimakofF expansion [51] on its 
fully polarized ground state. The quadratic Hamiltonian 
obtained (linearized in the fluctuations) is diagonalized 
to calculate the structure factor. The simulations were 
performed for systems of size L = 500,1000, 2000 with 
D = 3. 


The ground state of a fully F spin chain is a classical 
fully polarized state and the Hamiltonian governing sin¬ 
gle spin-flip excitations is given by a model of particles 
in random (Hartree-like) potentials and randomized tun¬ 
neling terms. Thus, one may expect these excitations to 
be localized. On the other hand, the failure of the RG 


















20 


10^ 

T . . . . . . . 

o 

o 

(a) IPR(a;) ~ td 

10^ 


10“^ 

- 

o 

o 





. 1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 

10 

“2 iQ-i 10° 1C 


10 “° 10 ^^ 10“2 10 ° 



2000 3000 



FIG. 16. (a) The finite-q structure factor Sq{u!) \q = 

(27r/25,7r)] as a function of frequency; (b) the fidelity / 
and (c) the noise frequency are plotted as a function of the 
RG steps for a two-dimensional strip of spins of lattice size 
500 X 6, with D — 3 a,t T = oo. There are ~ 50 clusters re¬ 
maining when these are of the size of the width of the spin net¬ 
work; these remaining clusters behave as a one-dimensional 
network of effective spins and produce the anomalous 
noise. Here we see that about 50 (red dashed lines) spins re¬ 
main when the maximum frequency a;~2xl0“^ and this 
approximately agrees with the frequency below which l/uj°‘ 
noise behavior is observed. Above this frequency, a ~ l/y^ 
tail is observed as in the case of the square network of spins. 




FIG. 17. (a) The inverse participation ratio is finite [seen by 
scaling with system sizes L = 500 (yellow), L — 1000 (blue) 
and L = 2000 (red) ] and scales with frequency as IPR(a;) ~ uj 
(dashed line fit); (b) spin-wave peak seen in the structure 
factor Sq{uj)-, (c) the peak occurs at frequency oJciq) ~ 
(red line fit corresponds to (d) the width (calculated 

as inverse height of the peak) scales as T{q) ~ q^ (red line fit 
corresponds to q^'^) 


suggests that the system is self-averaging and may pos¬ 
sess goldstone modes of a disorder-free ferromagnet with 
a spectral peak at frequency a;c(q) ~ It should be 
noted that the on-site potential and tunneling are exactly 
correlated in this system, thus, localization is not imme¬ 
diately guaranteed. The results in Fig. 17, 18 suggest 
that the system exhibits a bit of both; spin-waves exist 
with a central frequency Wc('Z) ~ but they are heav¬ 
ily damped with a decay rate F that scales in the same 
way, T{q) ~ q^. Thus, spin waves states are localized in 
the sense that the inverse participation ratio IPR [given 
by IPR = J tbe spin wave density 

at site i of the chain] is finite. If one interprets the fi¬ 
nite IPR as a consequence of an exponentially localized 
wave-function, then IPR ^ 1/^, where ^ is the localiza¬ 
tion length. Extracting the localization length from the 
IPR yields the scaling ^(w) ~ l/w; this is in agreement 


with the fact that the RG finds a dynamical exponent 
z = 1 in the limit of the fully F chain which yields a 
scaling between the cluster size n and the frequency uj as 
n(uj) ~ 

We note that this ‘localized’ spin-wave behavior is 
unique to one dimension. In two dimensions, an anal¬ 
ogous Holstein-Primakoff analysis reveals that the spin 
waves are sharp; the central frequency Wc q^ but the 
damping scales as F ^ g"* (see Fig. 19). Note that this 
is in agreement with the analysis of disordered ferro- 
magnets using an effective medium approximation dis¬ 
cussed in Refs. [41, 42] (the authors focus on the three- 
dimensional case, but their results are straightforwardly 
extended to two dimensions). The simulations for the 
two-dimensional case were performed for systems of lat¬ 
tice sites 20 X 20, 30 x 30, 40 x 40 and with disorder 
strength D = 3. 
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FIG. 18. An example localized spin-wave wave-function for 
the disordered purely F chain of size L = 500. 



FIG. 19. (a) The inverse participation ratio decreases to zero 
in the thermodynamic limit [seen by scaling with system sizes 
20 X 20 (yellow), 30 x 30 (blue) and 40 x 40 (red) ] ; (b) spin- 
wave peak seen in the structure factor Sq{u})', (c) the peak 
occurs at frequency tUc((/) ~ (red line fit corresponds to 
(d) the width (calculated as inverse height of the peak) 
scales as r(g) ~ q* (red line fit corresponds to g®'®®). 



FIG. 20. Infinite-temperature S^{uj) for classical spins of 
magnitude = 1/2 is plotted (blue) and compared with 
the quantum result (red) at infinite temperature. L = 15000, 
D = 3, gi = 0. 



FIG. 21. Comparison of the structure factor obtained for clas¬ 
sical spin chains using direct integration of Bloch-equations 
and RSRG rules. The direct simulation is computationally 
hard to perform over very long time-scales, but the agree¬ 
ment is very reasonable for the simulated frequencies. 


Appendix L: Classical simulations for T = oo 


The RG rules were derived using classical reasoning 
and can be applied to classical vector spins. Fig. 20 
plots (blue) the structure factor at infinite-temperature 


obtained from such a calculation with classical Heisen¬ 


berg spins of vector magnitude 


S 


1/2. The structure 


factor is very similar to the quantum result (red). We 
conclude that the anomalous dynamical exponent at in¬ 
finite temperature, that gives rise to anomalous diffusion 
is not a purely quantum phenomenon. 
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Appendix M: Comparison of classical RG with 
direct dynamical simulation 

The classical spin chain (size L = 1000) with ini¬ 
tial coupling distribution Po(|>^|) ~ |J| € 


[10“^, 10“^e®] is simulated using the classical Bloch- 
equations to obtain the structure factor. The result 
is compared with the classical-RG simulation result in 
Fig. 21. The simulations are for T = oo (the initial spin 
configuration is chosen entirely randomly). 


